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Deuteron life-time in hot and dense nuclear matter near equilibrium 
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We consider deuteron formation in hot and dense nuclear matter close to equilibrium and eval¬ 
uate the life-time of the deuteron fluctuations within the linear response theory. To this end we 
derive a generalized linear Boltzmann equation where the collision integral is related to equilibrium 
correlation functions. In this framework we then utilize finite temperature Green functions to eval¬ 
uate the collision integrals. The elementary reaction cross section is evaluated within the Faddeev 
approach that is suitably modified to reflect the properties of the surrounding hot and dense matter. 
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I. INTRODUCTION 


The complicated dynamics of heavy-ion collisions provides a great challenge for many-particle theory. In particular 
at intermediate energies where the elementary ingredients are rather well known - in terms of constituents and their 
respective interactions - the main problem arises from a sufficient description of the many-particle aspect. To provide 
single-particle distribution functions such reactions can be simulated on the basis of kinetic equations as, e.g., supplied 
by the Boltzmann-Uhlenbeck-Uhling (BUU) approach (see, e.g.. Refs. [|i|-^). 

However, the formation of light clusters such as deuteron^ helium, a-particles, etc., is an important phenomenon of 
heavy-ion collisions at intermediate energies, see, e.g.. Ref. 0]. Empirical evidence, including recent experimental data 
on cluster formation , indicate that a large fraction of deuterons can be formed in heavy-ion collisions of energies 
below E/A < 200 MeV. Also, during the expansion of the system the density can drop below the Mott-density of 
deuteron dissociation pT|-p^. 

The description of the formation of such bound states (clusters) during the expansion of hot and dense matter 
is not as well elaborated as the single-particle distribution. The main obstacle is that the formation of bound 
states requires the notion of few-body reactions within the medium. Even the simplest case, i.e., the abundances of 
deuterons that are determined by the deuteron formation via NNN dN {N nucleon, d deuteron) and break-up, 
dN —f NNN, reactions, requires a proper treatment of the effective three-body problem. Previous studies of the 
kinetics of deuteron production have utilized the impulse approximation to calculate the reaction cross section at 
energies above 200 MeV [^. For lower energies, viz. E/A < 200 MeV, the impulse approximation fails badly and a 
full three-body treatment of the scattering problem is necessary [Q. Furthermore, a consistent treatment of cluster 
formation in expanding hot and dense matter requires the inclusion of medium effects into the respective elementary 
reaction cross sections as has been done in the nucleon nucleon (NN)-case and proven to be substantial in BUU- 
simulations of heavy-ion reactions Q . Therefore we present an exact treatment of the three-body problem including 
medium modifications in mean-held approximation. 

The cluster formation during the expansion of hot and dense matter is driven by the collision term in the generalized 
linear Boltzmann equation (see, e.g., Ref.s |T1,@). Here we consider the huctuations of the deuteron distributions in 
hot and dense nuclear matter in a near-equilibrium situation. One important question in this context is the time scale 
of formation and disintegration processes that govern the evolution to chemical equilibrium. Within linear response 
theory we relate the reaction rates to the equilibrium correlation functions. To calculate the response coefficients it 
is then possible to apply the method of hnite temperature Green functions 

The essentials of the three-body problem for the isolated system are well known, see, e.g.. Ref. [|^. In the 
following we utilize the AGS-formalism suitably modified to treat the three-body problem within nuclear matter. 
To derive the proper AGS-type equations we use the self-consistent random phase approximation |22| extended to 
finite temperatures. For a numerical solution we rely on a separable representation of the NN-potential. This 
choice simplifies the problem considerably. A systematic investigation of separable parameterizations of “realistic” 
potentials has been pursued, e.g., in Refs. [^. We note that solutions of the three-body problem using “realistic” 
NN-potentials have been achieved, e.g., by the Bochum group 241, and the Bonn group in the framework of the 
IV-matrix approach 

In the following section we present the formalism to treat cluster formation in a linear approximation of the 
generalized Boltzmann equation. In Sec. HI we introduce the finite temperature Green function and derive a Faddeev- 
type equation that includes medium modihcations due to Pauli-blocking and energy-shifts. We relate the “in-medium” 
cross section to the collision term in the Boltzmann equation. Our numerical results are presented in Sec. IV and we 
summarize and conclude in Sec. V. 
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II. QUANTUM KINETICS AND BOUND STATE FORMATION 


The Hamiltonian of the Fermi system in question is given in terms of creation and annihilation operators, 

i? = + ^ E ^2(12,l'2')aI4a2mn, (1) 

11 / 121 ' 2 ' 

where a\ a satisfy the well known commutation relations. The indices 1,2,... collectively denote the quantum numbers 
(e.g., momentum, spin, isospin,...) of the particles 1,2,... The observed physical quantities will be expressed in terms 
of reduced n-particle occupation matrices (see, e.g.. Refs. 

/„(1... n, 1'... n'; t) = {a\, .. .a\,ai... a„)* = Tr{p(t) a\, ... a\,ai ... o„}, (2) 

where p{t) denotes the density matrix of the many-particle system. In case of equilibrium {p(t) = po) we use the 
notation /°(. ..) = (.. .)o. 

The reduced density matrices given above are particularly suited in the framework of a cluster decomposition . 
If clusters are treated in mean-field approximation we may introduce cluster wave functions pv and introduce bosonic 
(two-particle) operators that are given through 

^1/ = ^ aia2 (3) 

12 

and the h.c. . For the two-particle system of interest here pi, is given by the solution of the respective two particle 
Bethe-Salpeter equation with the eigenvalues (of bound or scattering states) ||^. This way it is possible to write, 
e.g., /2 in a cluster representation, viz. 


f2{v:v'\t) = {hl,KY. ( 4 ) 

For nuclear matter the conditions in the final stage of a heavy-ion collision may be such that formation of bound 
states is possible. This is indeed the case, when the density of the system is below the Mott density, e.g., of the 
deuterons, and formation will occur [pT[ - |T^ . Following a general density matrix approach as given in Ref. the time 
evolution of the distribution functions for nucleons fi{t) = with momentum p and deuterons f 2 {t) = fd{P,t) 

with momentum P reads for homogeneous matter 

dtfN{p,t) = Tr{p{t)i[H,nNp]} = -VN{p,t) +lNip,t), (5) 

dtfdiP.t) = TT{pit)i[H,ndp]} = -Vd{P,t) +Xd(P,t), (6) 

where njvp = oYfpdNp and ridp = b^dP^dP- The Vlasov terms P(t) describe the reversible time evolution and are 
related to time dependent Hartree-Fock calculations as, e.g., explained in Ref. p^ . 

The collision terms I{t) = +T^(t) correspond to the irreversible behavior and describe elastic scattering (E) 

and inelastic (reaction) processes (R), respectively, between the constituents of the system. The elastic processes do 
not change the internal quantum numbers of the particles. They determine the time scales for thermal relaxation. 
However, the inelastic processes that are related to excitation as well as to bound state formation and disintegration 
change the abundances of the components characterized by the internal quantum numbers and determine the time 
scale of chemical equilibration. The reaction relevant for the energy domain considered are due to photodisintegration 
I-yNN,dip,t) and nucleon deuteron break-up lNNN,Ndip,t) (and the reversed ones), i.e. 

P^{p,t) = I-iNN,d{p, t) + lNNN,NdiP, t) + . .. ( 7 ) 

Td'(-P) t) = Id,-yNNiP, t) + IdN,NNNiP, t) + . . . (8) 

Further reaction channels represented by the dots are given in Ref. p5| ]. Presently, we consider the three-particle 
processes. 

The collision integral that involve three-particle processes have been given in Born approximation or evaluated 
in impulse approximation fl^ . In both cases the influence of the surrounding medium on the elementary cross section 
that enter into the collision integrals has been neglected. This might not be sufficient for intermediate energy heavy- 
ion reactions as has been shown, e.g., for the NN collision rate in a BUU calculation of La on La |^. The approach 
presented here naturally respects medium modihcation in the break-up cross sections. Furthermore, in view of the 
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rather moderate energies reached, the effective three-body problem arising in this context is treated exactly in terms 
of properly generalized Faddeev-type equations. 

To be more specific we consider the situation where the collision rate is sufficiently high compared to the reaction 
rate, so that each components are close to their thermal equilibrium distributions. The small deviations of the chemical 
composition from equilibrium are then treated within the linear response theory. The time scale of the relaxation to 
chemical equilibrium is set by the reaction processes that will be considered in the following. In this case 

dtfNiP^t) = ( 9 ) 

dtff{P,t)=lf{P,t), (10) 

where we have introduced I^jp, t) = lNNN,Nd{p,t) and lf^{P,t) = IdN,NNN{P,t) for brevity. The collision terms 
on the rhs. of Eqs. (|^) and (|^) each contain gain and loss terms due to deuteron break-up or formation reactions. 
Collisions of higher clusters (e.g., dd) that require a suitable treatment of the effective four-body problem are left for 
further investigations. To evaluate the integral I^{P,t) we use linear response theory (see appendix A), 

= - X! {'ndp;ndP’){ndP>;ndP")~^ Sfd{P",t), (11) 

P'P" 


where Sfd{P, t) = fd{P, t) — fdiP) denotes the fluctuations from the equil ibrium distribution, and Udp = {b\pbdp)- The 
Kubo scalar product (A; B) that ap pears in Eq. dO) is given in Eq. (|All[) and its Laplace transform, i.e. the correlation 
functions (A (?7 —s- O’*");!?), in Eq. ( Al2| ). Following standard many-body techniques (see, e.g.. Refs. ]l^ , p^ , p^ ) the 
correlation function is evaluated using Green functions. 


{A{rjy,B)=- 


1 


duj 


1 


P J 2tt rj + iuj uj 


[G^B(a; -|- fO’’") — Gab{<^ — fO"*")] 


( 12 ) 


where Gab{z) is the analytic continuation of the Matsubara Green function GAB{Zfi) that will be discussed in the 
next section. 

For homogeneous matter where {hdp;ndp>) and {ndP'\ndP") are diagonal in momenta P, the response equation 
& is given by 

dtffiPA) = {ndpiv 0 +);hdp)(ndp;ndp)"^ Sf^{P,t) = — Sff{P,t)- ( 13 ) 

TdP 

The limit 77 ^ O’*" implied through Eq. (H) has to be taken after the thermodynamic limit. Here we have introduced 
the momentum dependent life-time (formation-time) of the deuteron fluctuations r^p in the surrounding medium, 
which is of central interest. Note that the disintegration (formation) of deuterons requires the explicit treatment of 
three-particle equations in nuclear matter, which will be derived in the following section. 


III. FINITE TEMPERATURE GREEN FUNCTION AND THREE-BODY EQUATIONS 


In order to evaluate Eq. (13) by use of Eq. we need to define the finite temperature Green function. To consider 
n particles embedded in a medium the n-particle Green function ^„(1... n, 1'... n') for equilibrium is defined by 

5^*'(1... n, F ... n') = -7 (T A„(t)Ai,(t))o, (14) 

where T implies Wick time ordering. The operators A„(t) are n-particle Operators, i.e. An(t) = 
ai(ti)... an{tn)\ti=...=tn=t taken at equal times, and in the Heisenberg picture 

A{t) = exp(iiLt) A exp{—iHt). (15) 

The Green functions given in Eq. (|lj) satisfy a hierarchy of equations, given, e.g., in Refs. 00- To arrive at 
equations that are solvable in practice for the n-particle problem, one has to truncate the hierarchy, which is usually 
done by introducing suitable approximations for the {n + l)st-particle Green function. For the three-particle problem 
this has been done, e.g., in Ref. in special cases. Within the self-consistent random phase approximation it is 
possible to arrive at equations that are already decoupled. This methods has been used for zero temperatures in 
Ref. 1221 and will be extended to finite temperatures here. To simplify the notation we use a matrix form in the 
following 
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( 16 ) 


The time evolution of the n-particle Green function is governed by a Dyson equation 


idtQ 


t-t' 

n 


+ J dtMi ‘ Gl, 


(17) 


The mass matrix * introduced in the above equation is given by 


M 


t-t 

n 


S{t — t)MnQ + M 


t-t 

n,irr. 


with 


(■^n,o-^n) (l--.n, l'...n') =Tr{po [An,H],Al, 

{l...n,l'...n) = -iTT{poT[A^,H]t[H,Al]i]}„r., 


(18) 


(19) 

( 20 ) 


where the index irr. indicates that all reducible parts should be omitted, where the index ± refers to odd (+) or 
even (—) numbern of fermions. The first term refers to an n-body cluster mean-held contribution fisl , ^^ whereas 
the second term is of dynamical origin and contains retardation. Up to the correlations of interest instantaneous and 
dynamical contributions are separated. The normalization is given by 

A 7 „(l...n,T...n') = Tr{po K,A|^,]±}. (21) 


In cluster mean-held approximation the term M.n,irr. will be neglected. In the Matsubara-Fourier representation the 
Green function is given by 


= —pT. 


( 22 ) 


For a fermionic system considered here is a fermionic or bosonic Matsubara frequency, depending on whether n is 
odd or even, resp. to preserve the Kubo-Martin-Schwinger boundary condition |17|-|l9|. 

Taking in Eq. (0) for As = 010203 and evaluating Eq. ( p^ ) in the independent particle approximation leads to the 
following Bethe-Salpeter equation for the three-particle Green function at hnite temperatures and densities, 

Gsiz^) = G^°\z^)+Rl^\z^)VsGsiz^), (23) 

which is the central input to derive Faddeev type equations in a medium. The notation will be explained in the 
following. The proper symmetrization is treated separately. The Green function of the noninteracting system is 

G^^\z^) = NsR^Hz^), (24) 


where 


(123,F2'3';z^) = 


^11'^22'^33' 

— £l — ^2 — £3 ’ 


-^ 3 ( 123 , 1'2'3') = (5i1'(522'<533'(/i/2/3 + / 1 / 2 / 3 ) 

= <5ii'i522'<533'(1 — fi — /j)(l — fk + gi^i + £j))- 


(25) 

(26) 
(27) 


Note that Eq. (27) is identical to Eq. ( p6[ ) for all permutations of ijk = 123. We use / = ! — /, and the Fermi 
one-particle function /(ei) = (exp[/3(ei — fi)] + 1)“^ and the Bose function g{uj) = (exp[/3((^ — 2/r)] — 1)“^ for the 
two fermion system. Here (3 is the inverse temperature of the system and /i is the chemical potential. In mean held 
approximation the single quasi-particle energy ei is given by 


£1 = 


3. 

2mi 


^HF 


(1), 


^HF 


(1) = ^[1/2(12,12)-U2(12, 21)] /2. 
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(28) 

(29) 







( 30 ) 


Note that = 0. The interaction kernel in V 3 in Eq. is given by 

3 

E3(123, T2'3') = E3^''^(123, T2'3'), 

fe=i 

F3('=Hi23, T2'3') = (!-/,- fj)V2{ij,i’f)Skk', (31) 

with ijk = 123 cyclic, and V 3 If we introduce a potential V 3 = we may instead of Eq. ( p^ write 

G^iz^) = Gf (z^) + Gf (z^)E3G3(^^), (32) 

which looks formally as the equation for the isolated case [p^. Using Eq. (|2^) we may write the potential V 3 = ^ 3 ^^^ 

in terms of (e.g., k = 1 ) 

U3(^^(123,1'2'3') = (l-/i+ 5(£2 + e3))"V2(23,2'3')5ii', (33) 

In Eq. (|^) we have already introduced the channel notation that is convenient to treat systems with more than 
two particles [ p 0 [ . 

If the pair and the odd particle are uncorrelated in channel 7, we may define a channel Green function G 3 '’^^(z^). 
In this case only the interaction within the pair of channel 7 is taken into account, viz. 

Y *'^2(wa) Gi(z^ - wa). (34) 


The summation is done over the Bosonic Matsubara frequencies uj\, A even, lo\ = TrX/{—i/3) + 2/r. The equation for 
the channel Green function is derived in the same way as for the total three particle Green function given in Eqs. (|2^) 
and (|3^). The result is 

g(^)(z^) = Gf (z^) + Gf (z^)U3(^)g(^)(z^). (35) 

Introducing the notation = V 3 — vY we arrive at the following equation for G 3 (z^) expressed through the 
channel Green functions GG\zfj), i.e. 

G3(z^) = g(^)(z^) + gG\z^)vYG3{z^). (36) 

Now we have set the necessary equations, i.e. Eqs. (p2|), (p5|) and (|^ ) to define a channel transition operator UajS 
fir finite temperature, 

G3 (z^) = 5^pGt\z^,) + G^^\z^)U^p{z^)Gf\z^). (37) 


In the zero density limit, this definition coincides with the usual definition of the transition operator with the correct 
reduction formula to calculate cross sections [p^. Inserting this definition into Eq. ( p^ and using Eq. (^) leads to 
an equation for the transition operator in medium, viz. 


= (1 - <5a/3)Gf^ vYg^^^U^p- 

i¥=oi 


(38) 


This is the AGS-type (or Faddeev-type) equation valid to treat three-particle correlations at finite temperatures and 
densities in mean-field approximation. Although this equation looks formally equal to that for the isolated system, we 
emphasize that vY as well as gG\ and hence is different from the isolated system due to the finite temperature 
and density of the surrounding matter, and therefore contains Pauli factors due to phase space occupation and self¬ 
energy shifts. This becomes transparent, if the definitions of the quantities appearing in this equation are inserted. 
Before doing that we define a transition channel operator , 


_ ^( 0 ) I ^(0) rpG) ^(0) 

'-'3 — '-'S +*-^3 23 '-^31 


inserting this equation into Eq. (Ba) leads to 


( 39 ) 
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( 40 ) 


and to This way it is possible to write a second, more useful version of the AGS-type equations 

Uo.0 = (1 - ^a/3)(A^34°V' + E (41) 


~f^a 


Here we have written the Pauli factors occuring due to the surrounding matter explicitly. Note, that through Eq. (H) 

is also medium dependent. To compare with our previous result Ill], we repeat the expressions for the low 

1/2 1/2 

density case. In this case, we may assume > 0, and therefore write an equation for = N^' U^pN^' 


Kp = (1 - <5„/5)i?f^ C/;^. 


~f^a 

The equation for the transition channel operator Tg = TVg^^Pg^^TVg^^ (and so for V3) is then 

rg*^'^) = Vg*(T') + 

Inserting all definitions the explicit form of the effective potential arising in this equation reads 
yg*(")(123, T2'3') = (1 - /i - -h + g{e, + £ 2 ))-'^' 


xy2(12, 1'2')(533'(I - /3 + 5(ei' + e2'))'/'(l - /i' - /2')'^' 
^ (1 - /i - / 2 )'/' 14 ( 12 , 1'20(1 - fv - / 2 ')'/^ 


(42) 


(43) 


(44) 

(45) 


where Eq. ( ^ ) holds for P f ■ Utilizing this approximation Eq. (^^ has been solved numerically using a separable 
ansatz for the strong nucleon-nucleon potential |14[| . 

We are now in the position to evaluate the correlation function ( 71(77 ^ O'*');?!) according to Eq. (12). The Green 
function needed in Eq. & is given by (see appendix B) 


^ E Tr {c/.o Gf (H^ + z.) U,, G(^p(z.)} 


(46) 




To perform the Matsubara summation that is present in Eq. ( ^^ , we now use the spectral representation of the Green 
functions that have been given for the quasi-particle approximation in Eqs. (|^) and (|^). The resulting expression 
for Gnj^upi^fj.) may be cast into the following form. 


g: 


(7) 


,(f^M)=4*E 


(47) 


(123|G-,o|c(^i]^)(c<^(]^|I/o^|123) 

123c + (^dP 4" ^c) ~ Eq 

X [/l/2/3/c5(i^^P^) - /l/2/3/c(l + ))] 

T(Hjy ^ ^ H^), 

with Eq = £1 -I- £2 + £ 3 , and c G {1'2'3'} depending on the channel 7 G {123}, resp. The terms in brackets of Eq. 
are usually referred to as Pauli factors of the gain and loss terms. Since presently we are interested in the time scale 
of fluctuations, we may consider, e.g., the loss term that is given by 


(^4p;^4p) = 4E l(123|Go^|(^i]),c)p/i/2/3/cff(4]{) 2^<5(4]{ +£e - Eo). 


(48) 


123c 


Eor identical particles that are considered here proper symmetrization has to be taken into account. To this end we 
connect the result given in Eq. ( |4^ ) with the transition matrix Tq, that is evaluated between properly symmetrized 
and normalized states (f>o, (t)d and again satisfies the equivalent three-body equation, see Refs. for the isolated 
case. It is given by 


(<).o|To|<^d) = y2E(123|Go^ 

W'dpP- 


(49) 
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To be more specific we also separate spin and momentum degrees of freedom and introduce amplitudes 


(mim2m3|A^o(kik2k3,kjvP;i^)|TOjTO2) = ((/)o|Td|</'<i)- 
Using this amplitudes the life-time may be written in the following way 


( 50 ) 


3! 


/ (fkid^k2cPk3 Tr(AfoPi7Wj) / 1 / 2 / 3 /E 2Tr6{E - Eq). 


(51) 


The total energy is given hy E = P^/2md + e + Ed and e = fc^/2m is the energy of the odd nucleon. The factor 1/3! 
prevents overcounting due to the six possible ways of arranging the identical particles among the three momenta, the 
trace is over spin projections only and pi is the initial spin density matrix. 

It is instructive to discuss several ways to recover the Born approximation and impulse approximations that have 
been used previously |]^,^. From Eq. ( ^ ) the lowest order iteration for the (on-shell) break-up amplitude U^q is 


C^70^E^3 -E^3 =^3 ’ 


(52) 


where the first term is referred to as impulse and the second as Born approximation. Replacing LUn by in 

E< 1 . 0) Ms to 


G 


(j)Bo 


‘(O ) = Tr( 




3 "-^3 

-u,), 


3,dP 


(53) 


A second possibility is to expand the Green functions into the spectral representation. In Eq. (^ ) the spectral 
representation of the Green functions leads to matrix elements (t/^p lUiyol'/o)- Then we may use the on-energy-shell 
relation 


(0o|c/o^l4p) = ^ 


which in turn after reinserting the spectral expansion leads again to Eq. (| 


(54) 


IV. RESULTS 


Although formally rather simple the integration over the momenta ki of Eq. (|51| ) is rather tedious, since momenta 
dependencies appear also in the Fermi functions /. In the low density approximation we may as indicated at the end 


of Sec. Ill, use the definition U = ^ which leads to ^Jvl *and 

-^3 — /l/2/3- 


(55) 


Hence the term / 1 / 2/3 in Eq. & will be absorbed in the redefinition of A4*. The additional Fermi functions 
(/i/ 2 / 3 )“^ appearing in the Nd channel due to the replacement A4o ^ Adg will be approximated in the following 
way (note p < /) 


fe/{flf2f3) — /e(1 + /l + /2 + /s) — fe 


(56) 


We may now introduce the in-medium break-up cross section cTg in the center of mass system, which coincides with 
the usual in-vacuum definition in the zero density limit 0- It is given by 


aPE) = 




|Vd - Vat I 


I 

3! 


d^p'd^' TT{MoP^Ml'^) 2TrS{E* - E*). 


(57) 


and \vd - Vat I is 
= I Am* + where 


The cross section is evaluated in the center of mass system introducing Jacobi coordinates p', q' 
the relative velocity of the incoming particles. The center of mass scattering energy is E* = 
we have used effective mass approximation for the nucleon self energy 1141. Due to the medium effects the deuteron 
binding energy changes, which is calculated consistently with the two-body input into the Faddeev equation that lead 
to the amplitudes A4 q. To evaluate the life-time of the deuteron in medium we introduce the cross section defined in 
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Eq. (|^) into Eq. The remaining integration is over the momentum k^v of the odd nucleon. The equation for 

the life-time of the deuteron is then given by 


‘dP 


4 



kN |vd - Vat I (To(^*) /(£*). 


(58) 


and £* = (3q/2 -|- P/2)^/2 to*. 

The cross section entering into Eq. is given in Fig. |^. We restrict the two-body channels to the dominant ones, 
i.e. ^5'o and ^Si — ^Di. For the separable ansatz we use the parameterization of Phillips The parameters are 
taken from Ref. |^, which lead to a good overall description of the elastic and break-up cross sections as well as 
the differential elastic cross section up to Eiab = 50 MeV |lJ,Q . To calculate the break-up cross section we use the 
optical theorem. 

The solid lines represent the isolated break-up cross section. As shown previously it reproduces the experimental 
data mil- The dashed lines show the break-up cross section, a^rp{Eiab) for densities n = 0.1, 1, 3, 5, 7 x 10 ^ 
fm“^, respectively, as a function of the laboratory energy Eiab- The Mott transition occurs at the density n ~ 8 x 10“^ 
fm“^. 

The medium effects significantly modify the isolated break-up cross section. Two qualitative features are observed. 
First, the break-up threshold is shifted towards lower scattering energies with increasing density of the nuclear matter. 
This kinematic effect is due to the decrease of the deuteron binding energy with increasing density (see Fig. ^). Second, 
the cross section increases considerably with increasing density. The maximum is enhanced by one order of magnitude 
for the largest density value considered. For densities larger than the Mott densities the deuteron disappears as a 
bound state. 

Also, it is instructive to see how the medium dependent cross section converges to the isolated one. At Eiab — 100 
MeV the deviation of the in-medium cross section from the free cross section in this model is in the order of 10%. 
From inspection of Figure ^ we conclude that the dominant changes in the cross section takes place at rather moderate 
energies, i.e. where the impulse approximation fails, and the Faddeev technique has to be used. 

The resulting deuteron life-time evaluated with Eq. ( |^ ) is shown in Fig. ^ The influence of the medium modification 
through the cross section is substantial, in particular for small deuteron momenta. For higher momenta the differences 
become smaller, as they should, since the medium effects on the single particle properties become smaller at higher 
momenta. Also, the difference between the Maxwell and the Fermi distributions are shown that are comparably small 
for the densities considered. 

In Fig. ||the width of the deuteron at P = 0 fm“^ is shown as a function of nuclear density. Respecting the medium 
effects in the cross section leads to a larger width of the deuteron of almost a factor of two near the Mott density. For 
comparison also the medium dependent deuteron binding energy is shown in the same scale. 


V. SUMMARY AND CONCLUSIONS 

As already expected from earlier results on the NN-cross section, we And that the deuteron break-up cross section 
{Nd NNN) is also substantially modified at finite densities and temperatures compared to the isolated one. The 
densities and the temperature chosen are expected to be typical values for the Anal stage of heavy ion collisions at in¬ 
termediate energies. To reach this conclusion we have extended the AGS-formalism to treat the effective three-nucleon 
problem in an environment of hot and dense nuclear matter. This has been achieved using the finite temperature 
Green-function method within the Dyson-approach. The three-body problem is then formulated in the cluster mean- 
field approximation, and the resulting AGS-equations are solved numerically for a separable NN-potential. For the 
isolated system the experimental data are reproduced within a few per cent. 

Within mean-field approximation the influence of the surrounding matter leads 1) to a shift of the self energy of 
the nucleon and deuterons and 2) to additional phase space factors due to Pauli blocking. This two effects are taken 
into account consistently as the three-body equations are solved. 

The influence of the medium on the break-up cross section is calculated in the center of mass frame for the three- 
particle system. It shows three important features: 1) The cross section increases with increasing densities, 2) the 
threshold energy shifts due to the decreasing binding energy of the deuteron at increasing densities, and 3) the effect 
of the medium becomes smaller at higher scattering energies. The influence of the medium is rather strong. Near the 
maximum and close to the Mott density the cross section increases almost one order of magnitude compared to the 
isolated one. We argue that this modification is also important in a complete treatment of the heavy ion reaction (as 
found for the NN-case Q). 

An important global quantity that governs the time scale of the deuteron formation is the life-time (i.e. width) of 
the fluctuations in the deuteron distributions in hot and dense nuclear matter. We have calculated the life-time of 



the deuteron fluctuations using either isolated or medium dependent cross section in the collision integral. We And 
that the life-time strongly depends on the type of cross section included in the evaluation. The difference between 
the use of the isolated cross section versus the medium dependent one amounts to almost a factor of three near the 
Mott density. This is in support of the statement that the medium modiflcations of the break-up cross section may 
lead to changes in the final outcome of the deuteron rate in heavy-ion collisions. Therefore the results presented here 
may be considered an important input for transport equations that are used to describe the dynamics in a heavy-ion 
collision as done, e.g., in Ref. Q. 

The formalism presented here is capable to be extended to effective n-particle equations and therefore to treat the 
formation of higher clusters than deuterons. In particular the formation of helium, triton, or/and a-particles is of 
special interest. In this context, the inclusion of the total momentum dependence of the in-medium cross section is 
necessary. 
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APPENDIX A: LINEAR RESPONSE AND THE COLLISION INTEGRAL 


In order to derive the exact relation for the collision integral with the full medium dependent cross section given in 
terms of properly defined medium dependent three-particle transition operators, we assume small fluctuations of the 
equilibrium distributions and utilize the linear response theory to treat the nonequilibrium aspect of the process. 

To do so consider first equilibrium. In this case the two-particle distribution function may be decomposed in the 
following way, reflecting the uncorrelated and the correlated part 32 , 


/°( 12 , T 2 ') = f°{l)f°{2)[SivS22' - Si2'S2v] + < 72 ( 12 , T 2 '), 

where in cluster mean-field approximation, 

52°(12, 1'2') = ^ ((12|<^,)(<^,|T2')-5.p(12)(5,p(1'2')U 9iE.), 


and d^p(12) = ^p.ki+k 2 '^(ki-k 2 )/ 2 ,i/i and ()ex denotes inclusion of exchange terms.Using Eqs. (M) and (| 


(Al) 

(A2) 
leads to 


/2 (u = g{E^)5^y , (A3) 

i.e., /2 ( 1 ^, is diagonal in the indices /i that will be used in the following. 

In the framework of linear response the one- and two-particle distributions may be characterized by the deviations 
from the respective equilibrium distributions via small fluctuations Sf, viz. 


/i(l;f)=/°(l)+<5/i(l;t), (A4) 

f2{v,v']t) = f2{v,v') + Sf2iiy,E;t). (A5) 

The relevant statistical operator properly including the one- and two-particle distributions is given by the gener¬ 
alized Gibbs state. 


Prei = ^ exp 


-f3{H - fiN) - - / 3 ^F 2 (i/, v\t)Sn 2 {v, v') 


(A 6 ) 


where the operators describing density fluctuations of the one- and two-particle distributions appearing in Eq. ( [A^ ) 
are defined by 


(5ni(l) = a\ai - (aiai)o = ni(l) - /i°(l), 

5n2{y, u') = bl,b„ - {blX)o = n 2 {v, v') - v'). 


The Lagrange parameters Fi (t) and F 2 (t) of Eq. (|A 6 [) are determined by the consistency relations 


(A7) 

(A 8 ) 
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/k(0 = Tr{pi.ei(t) n„} 


(A9) 


whe re k collectively denotes the quantum numbers of the one- or two-particle operators. Linearizing Eqs. (A4) and 
([A^), resp., with respect to leads to 


= X! (AlO) 

k' 

where we have introduced the Kubo scalar product 

{A]B) = ^ j dr Tr{po A(-iT) B}. (All) 

0 

Furtheron we use the Laplace transform 

{A{rj)-B)= f dte'^^ (A12) 

J —OO 

For the case of the deuteron density fluctuation, Fq. ( |A^ ) reads 

Sridp = h\pbdp — g{Edp), (A13) 

so that 

{SndPy=Sff{PA). (A14) 

For small fluctuations 6f^{P,t) the response parameters Fdp{t) are small so that after linearizing we obtain the 
explicit relation 

5ff{PA) = ^ f3FdP'(t) {SndP'iSudp). (A15) 

P' 

The nonequilibrium statistical operator has the form (see e.g. p^), 

pit) = Prelit) - lim f dt' Uit, t'){i[H, Prelit')]+dt'Prelit')} uit', t)- (A16) 

7)^0+ J 


The occupation of the bound states tends to reach the equilibrium value due to the reactions within the system. After 
linearization with respect to Fdpit) and neglect of the explicit time dependence (Markov limit) in the integral, viz. 
Prelit') — Prelit), we obtain 


dt SmP, t) = {i[H, ndp])^ = Ifpit), 

where we have introduced I^pit) = IdN,NNNiP,t) for brevity. Evaluating Eq. ( AlTf ) leads to 
lfiP,t) = Tr{p(t) i[H,ndp]} = ^f3FnP'it) [indp'^fidp) - {hdP'^ridp)] ■ 


(A17) 


(A18) 


Note , tha t for homogeneous matter [ridp' ,nnp\ = 0, and than iridp'ifidp) = 0, which finally leads to Eq. dO) by using 
Eq. ( [Al^ ). 


APPENDIX B: EVALUATION OF THE CORRELATION FUNCTION 

For the one-particle occupation rik = i [H, rik] we get 

fik = i (1|^) ^ 2 (^ 2 ,1'2') a\a 2 a 2 'ap + h.c. (Bl) 

121 ' 2 ' 


In the three-particle space this leads to 
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(B2) 


n^k'^ = ^ a\a\a\aya 2 ’av + h.c., 

1231'2'3' 

where we have introduced the third particle in the channel 7 , and replaced V 2 ( 12 ,1'2') (3|3') ^ t^*'^^(123,1'2'3') (for 
7 = 1), which has been defined in the previous section. Further we use V 2 ( 12 ,1'2') = —V 2 ( 21 ,1'2') etc. The projection 
operator is given by (e.g. 7 = 1 ) 


= |fc23)(fc23|. (B3) 

The proper symmetrization of the final result will be treated separately within the framework of the three-body 
formalism. 

For n^p = ^dpb\p in channel 7 the time derivative is given by 

nl^P =—i (123|A^pF3^^^|l'2'3') a{a 2 a 3 a 3 'a 2 'ai'+ h.c., (B4) 

1231'2'3' 

where we have introduced ( 7 = 1 ) 


A^ =2 |W)(l(^dp|. (B5) 

To keep the notation transparent, we introduce a particle-hole (p/i)-basis, i.e. |123) 0 |I23) and extent the potential 
in the following way, 




'3, 


f(7) _ I 

^ 6,a - 'a 




(B 6 ) 


with a = {fc, dP}. 

Due to the operators occuring in Eq. (B4) the correlation function {n;n) is related to the 12-point Green function 
G 6 (D^) via Eq. dl)- Using Eqs. (^ ) and ( p4|) the Green function Gn^uki^ti) reads 


= Tr{r (]2 Gein^)Ilpk + (D. ^ -D.), 


(B7) 


and Ilph exchanges all particle with hole indices. The contribution with is due to the (h.c.) that appear 

in Eq. (B4). If only three-particle correlations are considered the full 12-point Green function is given by 


Ge{%) = -i- ^G3(D^ + za) G G3(za), 


(B 8 ) 


where explicit Matsubara summation has been introduced. Then summation over the block indices can be performed 
so that the resulting traces are in three-particle space only. Eq. (B7) simplifies to 


E Tr { G3(D^ + za) G3(za)} 


-l-(Dy —Dj,). 


(B9) 


In Eq. (|B9[) the term Ai"’'^G 3 (zA)Aa^'' appears. For a = k and e.g. 7=1 this is given by 

(I23|A[,^^G3A[,^P'2'3') = G3(/c23,/fc2'3'). 


. (7) 


(BIO) 


For a = dP we first consider using Eq. (|^, i.e. 


G3A(')) = G^^'>K^2p + 


dP' 


(BII) 


To evaluate this expression we introduce the spectral decomposition of Gg 


(7) 


Gi’FA) = Z 


iZ’xZ’i 


(ZA - £ 7 ) - E, 


(7) 


(l-/7+5(Pi^))). 


(BI2) 
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The wave function given in Eq. (B12) is a direct product of noninteracting wave functions of the pair and the odd 
particle, e.g. = \lLpn)- Through orthogonality, only the bound state part with momentum P contributes. 

Similar arguments hold for the projection from the left side. We will denote this cluster Green function by G^^dP- 
Note that through Eq. (|B11) etc. the full scattering solution is taken into account. Eq. (B9) then reads 


^AdpAdp 


+ (Ojy <—> 


v; 


( 7 ) 


3 1^3,0+ ^ a ) V3 

-n^). 


(7) 


G: 


i,dp(2:A)| 


(B13) 


Since we are only interested in the break-up reaction we have introduced an extended notation also for the second 
Green function appearing in Eq. ( Blj ), G 3 0 denotes the full Green function that describes the break-up situation. 
The related diagram is given in Fig. |^. The Born approximation is given by Fig. ^ 

It is important to note, that through Eq. ( 0 ) we are only interested in the discontinuity of on the real 

axis, which leads to energy conservation. Therefore we eventually need to consider only the on-energy-shell limit. It is 
now possible to introduce the break-up operator C /07 given in the Section III. Using the on-energy-shell requirements. 


G 07 GI = 


G^oGf = Ug'^G 


3,05 


(BI4) 

(BI5) 


these equations lead directly to Eq. ( |46|) and Eq. ([47|). For small deviations from equilibrium we may consider, e.g., 
the contribution of the loss term (first term in brackets of Eq. to determine the scale of the life time of the 

fluctuations. The resulting expression for the deuteron distribution in hot and dense nuclear matter is then given by 


'''dP = i^ndp'jSfidp) g (Edp). 


Here we have used that in ladder t-matrix approximation 


iridp^ndp) = g{Edp). 


(BI6) 


(BI7) 


This completes the deviation of Eq. 
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FIG. 1. Break-up cross section at temperature T = 10 MeV. Free cross section is shown as solid line and reproduces the 
experimental data, see Ref. Qj. Other lines are due to different nuclear densities, see text. 
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FIG. 2. Width of the deuteron in hot and dense nuclear matter at T = 10 MeV depending on the nuclear density and P = 0 
fm“^. The full triangles show the full calculation, the empty triangles show the one that uses the vaccuum cross section for 
the break-up reaction, free masses and deuteron binding energy. The deuteron binding energy is shown by the full diamonds. 
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FIG. 3. Momentum dependent life-time Tdp of the deuteron in hot and dense nuclear matter at T = 10 MeV. Upper two 
curves are without medium modified cross section using the solid line of Fig. |^, lower two lines with medium modifications 
for comparison. Solid and dashed-dotted lines uses Fermi distributions, the long dashed and short-dashed lines use classical 
distribution functions. 
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FIG. 4. Pictorial demonstration of the Born approximation for The dots indicate the Potential V 2 . Exchange and 

rearrangement channels have to be added for a full treatment. 



FIG. 5. Approximation of Ge with full treatment of the intermediate three-particle Green function. 
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